(12) INTERNATIONAL APHBTATION PUBLISHED UNDER THE PATENT 



(19) World Intellectual Property 
Organization 
International Bureau 





(43) International Publication Date 
11 March 2004 (11.03.2004) 



PCT 



(10) International Publication Number 

WO 2004/021028 Al 



(51) International Patent Classification 7 : G01R 33/563, 
A61B 5/026, 5/055 

(21) International Application Number: 

PCT/US2003/026580 

(22) International Filing Date: 26 August 2003 (26.08.2003) 

(25) Filing Language: English 

(26) Publication Language: English 

(30) Priority Data: 

60/406,040 27 August 2002 (27.08.2002) US 

(71) Applicant (for all designated States except US): 
KENNEDY KRIEGER INSTITUTE [US/US]; 707 
North Broadway, Baltimore, MD 21205 (US). 

(71) Applicants and 

(72) Inventors: VAN ZUL, Peter, C, M. [NL/US]; 7712 Mill- 
stone Ct, Ellicott City, MD 21043 (US). LU, Hanzhang, 
(nmi) [CN/US]; 4403 Falls Bridge Dr., Apt. F, Baltimore, 
MD 21211 (US). GOLAY, Xavier, (nmi) [CH/US]; 1129 
William Street, Baltimore, MD 21230 (US). 



(74) Agent: KOCOVSKY, Thomas, E.; Fay, Sharpe, Fagan, 
Minnich & McKee, LLP, 1100 Superior Avenue, Seventh 
Floor, Cleveland, OH 44114-2579 (US). 

(81) Designated States (national): AE, AG, AL, AM, AT, AU, 
AZ, BA, BB, BG, BR, BY, BZ, CA, CH, CN, CO, CR, CU, 
CZ, DE, DK, DM, DZ, EC, EE, ES, FI, GB, GD, GE, GH, 
GM, HR, HU, ID, IL, IN, IS, JP, KE, KG, KP, KR, KZ, LC, 
LK, LR, LS, LT, LU, LV, MA, MD, MG, MK, MN, MW, 
MX, MZ, NI, NO, NZ, OM, PG, PH, PL, PT, RO, RU, SC, 
SD, SE, SG, SK, SL, SY, TJ, TM, TN, TR, TT, TZ, UA, 
UG, US, UZ, VC, VN, YU, ZA, ZM, ZW. 

(84) Designated States (regional): ARIPO patent (GH, GM, 
KE, LS, MW, MZ, SD, SL, SZ, TZ, UG, ZM, ZW), 
Eurasian patent (AM, AZ, BY, KG, KZ, MD, RU, TJ, TM), 
European patent (AT, BE, BG, CH, CY, CZ, DE, DK, EE, 
ES, FI, FR, GB, GR, HU, IE, IT, LU, MC, NL, PT, RO, 
SE, SI, SK, TR), OAPI patent (BF, BJ, CF, CG, CI, CM, 
GA, GN, GQ, GW, ML, MR, NE, SN, TD, TG). 

Declaration under Rule 4.17: 

— of inventorship (Rule 4. 1 7( iv)) for US only 

[ Continued on next page J 



(54) Title: MICROVASCULAR BLOOD VOLUME MAGNETIC RESONANCE IMAGING 




k-spc_mem 
40 



Blood Tl 
meas seq. 
62 



00 



10 




Recon__proc 




Img_mem 




Vid proc 


44 




46 




50 



f^l (57) Abstract: A magnetic resonance imaging system includes a magnetic resonance imaging scanner (10) that performs an in- 
version recovery magnetic resonance excitation sequence (70) having a blood-nulling inversion time (60) determined based on a 
^ blood TI value appropriate for a selected magnetic field and blood hematocrit, whereby magnetic resonance of blood is substantially 
O nulled. The inversion recovery excitation sequence (70) includes an inversion radio frequency pulse (74) applied with a small or zero 
slice- selective magnetic field gradient pulse to avoid inflow effects, and an excitation radio frequency pulse (80). The inversion pulse 
(74) and excitation pulse (80) are separated by the inversion time (60). The magnetic resonance imaging scanner (10) subsequently 
performs a readout magnetic resonance sequence (72) or spectroscopy sequence to acquire a magnetic resonance signal from tissue 
^ other than the nulled blood. A reconstruction processor (44) generates a reconstructed image from the acquired magnetic resonance 
^ signal. 
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MICROVASCULAR BLOOD VOLUME MAGNETIC RESONANCE IMAGING 



Background of the Invention 

The following relates to the diagnostic imaging arts. It finds particular 
application in non-invasive measurement by magnetic resonance imaging of cerebral blood 
volumes, and will be described with particular reference thereto. However, it also finds 
5 application in measurement by magnetic resonance imaging of blood volumes in other tissues. 

In-situ measurement of blood volume is useful in various clinical, diagnostic, 
and research applications. Local cerebral blood volume changes, for example, correlate with 
local neuronal activity in the brain. Cerebral blood volume measurements during physiological 
stimulation thus provides a tool for functional studies of brain activity. Cerebral blood volume 

10 measurements can also provide information about impaired and/or damaged tissue in stroke 
victims, as well as about lesions in many disorders, including, but not limited to cancer, 
vascular disorders, and the like. Blood volume imaging of other organs besides the brain can 
similarly provide functional and diagnostic data that is useful in clinical studies, diagnoses, and 
tests (e.g. stress tests or tests of vascular compliance). 

15 Various imaging modalities have been used to measure blood volume 

characteristics, including positron emission tomography, single photon emission computed 
tomography, and magnetic resonance imaging. Most of these are invasive approaches, in which 
the subject is administered a suitable contrast agent that selectively enhances blood contrast in 
the selected imaging modality. For magnetic resonance imaging, various paramagnetic contrast 

2 0 agents are commonly used for this purpose. The requirement of an administered contrast agent 
is a substantial disadvantage of these techniques. 

Magnetic resonance imaging of blood oxygenation level dependence (BOLD) 
is a non-invasive technique for indirectly measuring blood volume. In this technique, blood 
hemoglobin is used as an endogenous contrast agent. In one BOLD imaging approach, 

2 5 magnetic resonance imaging is performed as a function of physiological stimulation that causes 
changes in blood oxygenation level. Blood volume is estimated from BOLD measurements by 
making assumptions pertaining to other parameters that affect blood oxygenation level, such as 
blood flow. Hence, BOLD does not provide a direct measure of the blood volume. 

A disadvantage of both the invasive techniques and the BOLD techniques as 

30 applied to blood volume measurement is that these existing techniques generally do not 
differentiate between blood in large blood vessels, on the one hand, and perfused blood in 
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small capillaries or other microvessels, on the other hand. The blood volume in larger blood 
vessels is principally controlled by sympathetic regulation. In contrast, blood volume in 
microvessels having typical diameters of less than about 200 microns tends to vary to maintain 
local homeostasis or in response to chemicals such as vasodilators or vasorestrictive 
5 compounds. Consequently, the blood volume of microvessels responds to physiological 
perturbations such as local neuronal activity. For functional magnetic resonance imaging, the 
volume of blood in the microvessels is typically of principle interest, while the blood signal 
from larger blood vessels is interfering and thus undesirable. On the other hand, total 
microvascular plus macrovascular blood volume may change in some diseases, including but 
1 0 not limited to for instance arteriovenous malformations. 

The present invention contemplates an improved apparatus and method that 
overcomes the aforementioned limitations and others. 

Brief Summary of the Invention 

According to one aspect, a magnetic resonance imaging method is provided. A 
15 blood-nulling magnetic resonance excitation sequence is performed that substantially nulls a 
magnetic resonance signal from blood. Subsequent to the performing of the blood-nulling 
magnetic resonance excitation sequence, a readout magnetic resonance sequence is performed 
to acquire a magnetic resonance signal from tissue other than the nulled blood. 

According to another aspect, a magnetic resonance system is disclosed. A 
2 0 blood nulling means is provided for performing a blood nulling magnetic resonance excitation 
sequence that substantially nulls a magnetic resonance signal from blood. A readout means is 
provided for performing a readout magnetic resonance sequence to acquire a magnetic 
resonance signal from tissue other than the nulled blood, the readout means operating 
subsequent to operation of the blood nulling means. 
25 One advantage resides in measuring the parenchymal vascular space 

occupancy, which is more sensitive to physiological perturbations than is the total vascular 
volume which includes the large blood vessels. 

Another advantage resides in providing measurements of the absolute blood 

volume. 

3 o Yet another advantage resides in providing images with substantially nulled 

blood magnetic resonance signal, i.e. a blood signal reduction sufficient to have MRI signal 
remaining that is predominantly from other tissues. 

Numerous additional advantages and benefits will become apparent to those of 
ordinary skill in the art upon reading the following detailed description of the preferred 

3 5 embodiments. 
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Brief Description of the Drawings 

The invention may take form in various components and arrangements of 
components, and in various process operations and arrangements of process operations. The 
drawings are only for the purpose of illustrating preferred embodiments and are not to be 
5 construed as limiting the invention. 

FIGURE 1 diagrammatically shows a magnetic resonance imaging system for 
measuring vascular space occupancy. 

FIGURE 2 diagrammatically shows a magnetic resonance imaging pulse 
sequence employing inversion recovery to null the blood signal. In FIGURE 2, the time 
1 0 intervals TI, TE, TR are not drawn to scale. 

FIGURE 3 shows a block diagram of a processor for computing a blood 
volume at rest and a blood volume change occurring responsive to a physiological perturbation. 

FIGURE 4 shows a block diagram of a processor for computing absolute blood 
volume using a reconstructed image acquired with the blood signal nulled by an inversion 
1 5 recovery preparation. 

FIGURE 5 shows a block diagram of a processor for computing a difference 
image having enhanced contrast for vascular space occupancy abnormalities. 



Detailed Description of the Preferred Embodiments 

With reference to FIGURE 1, a magnetic resonance imaging scanner 10 
includes main magnet coils 12, which are preferably superconducting coils, although resistive 
2 0 main magnet coils or a permanent magnet can also be employed. The main magnet coils 12 are 
energized to generate a substantially uniform main magnetic field in an examination region 14. 
Magnetic field gradient coils 16 produce gradients in selected spatial directions to spatially 
encode magnetic resonances that are generated by energizing a radio frequency coil 18. In 
FIGURE 1, a whole-body radio frequency coil 18 is shown; however, local coils such as head 

2 5 coils, phased radio frequency coil arrays, SENSE coils, and the like can be used instead of or in 

conjunction with the whole-body radio frequency coil 18 to excite magnetic resonances and/or 
to detect magnetic resonance echoes. 

A magnetic resonance sequence controller 30 coordinates and controls a radio 
frequency transmitter 34 that is coupled to the whole-body radio frequency coil 18 or another 

3 0 radio frequency coil to excite magnetic resonance echoes, and controls magnetic field gradient 

controllers 32 coupled to the gradient coils 16 to spatially encode the excited magnetic 
resonance echoes. One or more radio frequency receivers 36 coupled to the whole-body radio 
frequency coil 18 or another radio frequency coil detects, demodulates, and digitizes the 
magnetic resonance echoes and stores digital magnetic resonance samples in a k-space memory 
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40. A reconstruction processor 44 performs a Fourier transform-based image reconstruction or 
other type of image reconstruction to generate one or more reconstructed images from the 
stored k-space magnetic resonance samples. 

The reconstructed images are stored in an image memory 46, processed by a 
5 video processor 50 and displayed on a user interface 52, transmitted over a local computer 
network or the Internet, or otherwise processed. Preferably, the user interface 52 includes a 
display, printer, or other output device that allows a radiologist or other operator to view, 
render, or otherwise manipulate the reconstructed images. Moreover, the user interface 52 
preferably enables the radiologist or other operator to communicate with the magnetic 
10 resonance sequence controller 30 to create magnetic resonance imaging sequences, modify 
imaging sequences, execute imaging sequences, or otherwise control the magnetic resonance 
imaging scanner 10. 

With continuing reference to FIGURE 1, for clinical or diagnostic applications 
in which the blood volume is of interest, a measure of the vascular space occupancy (VASO) is 
15 performed. The approach uses an inversion recovery magnetic resonance excitation sequence 
having an inversion time 60 optimized for substantial nulling of the blood signal based on a 
blood Tl value particular to a certain magnetic field strength and hematocrit range. The 
inversion recovery magnetic resonance excitation sequence substantially nulls the signal from 
blood so that the magnetic resonance signal corresponds predominantly to tissue, with 
20 substantially negligible contribution of blood. This approach takes advantage of a substantial 
independence of the Tl value of blood on blood oxygenation, which allows the inversion 
recovery blood nulling to be operative for different types of vessels (arterial, arteriolar, 
capillary, venular, and venous) as well as throughout a physiological perturbation cycle. 

A suitable value for the inversion time (Tl) 60 for nulling blood can be 

2 5 obtained in a number of ways. In one suitable approach, a blood Tl measuring sequence 62 is 

applied by the sequence controller 30 to measure a Tl value of a representative blood sample 
using a blood perfusion apparatus, and the inversion time (Tl) 60 is computed therefrom. 
Alternatively, a major blood vessel can be identified in a reconstructed image, and a direct 
measurement of the Tl value of blood inside the identified blood vessel is obtained using 

3 0 conventional magnetic resonance imaging sequences. In another approach, as it has been found 

in practice that the Tl of blood is generally not subject-dependent for most human subjects 
having substantially normal blood, suitable inversion times (Tl) can be computed from 
measurements of a representative blood sample of substantially normal human blood or of 
animal blood of suitable mammalian origin and collected in a table 64, such as Table I and 
35 Table II contained herein. These values are generally suitable for substantially normal human 
blood, and include inversion time dependence upon the main magnetic field strength (Table I 
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providing values for a 1.5T field, and Table H providing values for a 3.0T field) and on the 
sequence repeat time (TR). Rather than employing tabulated data, the inversion time 60 can be 
related to magnetic field strength, repeat time TR, and optionally other parameters by an 
empirical functional relationship or other suitable relationship. The tabulated TI values are 
5 preferably used as a guideline for determining an inversion time TT that substantially nulls 
blood so that signals are predominantly from tissue. The tabulated values do not exclude other 
TI values that can also accomplish blood nulling. 



10 



TR (ms) 


optimal TI (ms) 


remaining gray signal 


1000 


409.4499457 


0.039848653 


1500 


551.5765181 


0.07104825 


2000 


659.2454703 


0.100852361 


2500 


738.9402994 


0.126845436 


3000 


796.8489702 


0.148292976 


3500 


838.3345867 


0.165337186 


4000 


867.7421743 


0.178518564 


4500 


888.427932 


0.188505318 


5000 


902.8982691 


0.195951912 


5500 


912.9810254 


0.20143462 


6000 


919.987164 


0.205430439 


6500 


924.8460607 


0.208318508 


7000 


928.2112602 


0.21039169 


7500 


930.5397637 


0.21187148 


8000 


932.1498927 


0.212922716 


8500 


933.262774 


0.213666531 


9000 


934.0317304 


0.214191053 


9500 


934.5629341 


0.214559874 



Table I - Optimal inversion times (TT) for blood at 1.5 Tesla magnetic field 
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20 
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TR (ms) 


optimal Tl (ms) 


remaining gray signal 


1000 


424.6207165 


0.029639854 


1500 


583.5056824 


0.054873106 


2000 


711.232168 


0.080671988 


2500 


812.0982526 


0.1047927 


3000 


890.5574158 


0.126155559 


3500 


950.8320255 


0.144394456 


4000 


996.6759107 


0.159557069 


4500 


1031.270597 


0.171909777 


5000 


1057.21771 


0.18181485 


5500 


1076.588295 


0.189656765 


6000 


1090.998208 


0.195801064 


6500 


1101.689382 


0.20057404 


7000 


1109.605744 


0.204255213 


7500 


1115.458796 


0.207077192 


8000 


1119.781556 


0.209229423 


8500 


1122.971521 


0.210863669 


9000 


1125.324122 


0.212099926 


9500 


1127.058395 


0.21303208 



Table II - Optimal inversion times (Tl) for blood at 3.0 Tesla magnetic field 

An inversion time adjustment is optionally performed to compensate for 
deviations in blood Tl value resulting from abnormal hematocrit values, sickle-cell 
pathologies, or another blood abnormality of a specific imaging subject. In one suitable 
adjustment technique, the blood-nulling inversion recovery imaging sequence is performed for 
several inversion times around the blood nulling inversion time selected from Tables I and II, 
and the blood nulling inversion time 60 is selected as the inversion time providing substantially 
negligible image signal from a large blood vessel A data flow line 66 in FIGURE 1 
corresponds to selecting the blood nulling inversion time 60 based on magnetic resonance 
measurements performed using the scanner 10. 

FIGURE 2 diagrammatically shows a suitable magnetic resonance imaging 
pulse sequence employing a blood-nulling inversion recovery magnetic resonance excitation 
sequence 70 to null the blood signal, and an exemplary single shot echo planar imaging readout 
72. Note that the echo planar imaging readout is exemplary only, and does not exclude 
employing additional or other magnetic resonance imaging, magnetic resonance spectroscopy 
or localized spectroscopy detection schemes. For example, the readout sequence 72 can be a 
single-shot imaging sequence, a single-shot echo planar sequence, a multi-shot imaging 
sequence, a spectroscopy sequence, a multiple slice image, a one-dimensional, two- 
dimensional, or three dimensional spatial encoding sequence, a fractional k-space acquisition 
sequence, a spin echo readout sequence, a gradient echo readout sequence, or the like. 



10 



15 



WO 2004/021028 ^^CT/US2003/026580 



7 

Loolcing first at the inversion recovery excitation sequence 70, an inversion 
pulse 74 is applied to invert the spins. Preferably, the inversion pulse 74 is a 180° pulse 
implementing a 180° flip angle for the spins. In this case, the inversion time TI 60 is 
substantially representative to induce proper blood signal reduction based on the Tl value of 
5 blood for the particular experimental conditions (field, hematocrit, oxygenation, etc). It is also 
contemplated, however, to use an inversion pulse having a flip angle greater than 90° but other 
than 180°, in which case the appropriate inversion time TI 60 is readily computed from the 
proper Tl and acquisition parameters such as the repetition time TR. In general, the inversion 
time TI 60 is selected as a time during which the longitudinal component of the flipped spins of 
10 blood decay from the flip angle to the crossover or null position. The null condition 
corresponds to a zero-crossing point of the longitudinal spin component; as the longitudinal 
spin component decays from the flipped or inverted alignment back toward the normal, 
non-inverted alignment it passes through a point where the longitudinal component passes 
substantially through zero, that is, the longitudinal spin component through a substantially zero 
15 crossing point. Substantially zero is understood to correspond to a substantially negligible 
blood signal such that the acquired magnetic resonance signal predominantly contains signals 
from tissues other than blood. 

Preferably, the inversion pulse 74 is not accompanied by a spatial encoding 
magnetic field gradient pulse or is accompanied by a relatively small spatial encoding magnetic 
20 field gradient pulse. This ensures that the spins of blood throughout the subject region of 
interest reach the null condition after the inversion time delay TI 60. In particular, the blood 
nulling is independent of blood flow since the blood-nulling inversion pulse 74 is spatially 
non-selective or selects a relatively large region. Thus, flowing blood that flows into the slice 
of interest at the time of excitation or at the time of readout is nulled appropriately. 
25 At the blood null condition, that is, after a delay time corresponding to the 

inversion time TI 60, an excitation pulse 80 is applied in conjunction with a slice-selective 
magnetic field gradient pulse 82 to excite spins in a selected slice of the subject region of 
interest. The excitation pulse 80 is preferably a 90° excitation pulse having a flip angle of 90°; 
however, an excitation pulse with other than a 90° flip angle is also contemplated. Because the 
3 0 blood is in a null condition at the time the excitation pulse 80 is applied, negligible magnetic 
resonance signal is excited in the nulled blood by the excitation pulse 80. Tissue such as fat, 
gray and white brain tissues, and the like generally have a different Tl value from that of 
blood, and so these tissues are not at a null condition at the time the excitation pulse 80 is 
applied. Hence, the excitation pulse 80 excites magnetic resonance predominantly in tissue. 
35 The exemplary single shot echo planar imaging readout 72 samples the 

magnetic resonance excited in the tissue by the excitation pulse 80. The illustrated single shot 
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echo planar imaging readout 72 is a conventional readout including a 180° spin refocusing 
radio frequency pulse 84 and slice-selective gradient pulse 86 that create a spin echo in a 
selected slice at a time-to-echo interval TE after the excitation pulse 80. A series of 
phase-encoding magnetic field gradient pulses 88 and a generally sinusoidal read magnetic 
5 field gradient waveform 90 step through a grid of k-space values in the selected slice while the 
radio frequency receiver 36 of FIGURE 1 performs sampling 92 of the spin echo. The magnetic 
resonance signal k-space samples are stored in the k-space memory 40 of FIGURE 1 and are 
processed to produce a reconstructed image representation as described previously. 

The single shot echo planar imaging readout 72 shown in FIGURE 2 is 

1 0 exemplary only. Substantially any type of magnetic resonance imaging or magnetic resonance 
spectroscopy readout or acquisition sequence can be employed. For imaging applications, the 
readout sequence may be a fast, single-shot sequence that acquires at least one slice per 
excitation. For example, the exemplary single shot echo planar imaging readout 72 can be 
replaced by a short gradient echo readout sequence. An echo time (TE) of the readout shown in 

15 FIGURE 2 is preferably kept short so that the readout is performed while the blood remains 
substantially in the nulled condition and to minimize contributions such as those from the 
blood-oxygen-level-dependent (BOLD) effect. Optionally, one or more additional inversion 
pulses are applied to keep the blood close to the null condition during longer readout 
sequences. 

2 0 Although the inversion recovery blood-nulling sequence 70 is preferred, other 

pulse sequences can be employed to substantially null the blood signal so that the magnetic 
resonance is predominantly due to tissues other than blood. For example, the use of a long echo 
time (TE) at high magnetic field produces a substantially blood-nulled signal; however, in 
sequences of this type the extravascular blood oxygenation level dependence (BOLD) 

2 5 contribution is large, complicating determination of other physiological parameters such as 

cerebral blood volumes (CBV) and the like. In one suitable approach for overcoming this large 
BOLD effect, the CBV contribution is quantified from interpolation of the TE dependence 
from a series of TE values at sufficiently long TE to avoid intravascular contributions that 
would occur at shorter TE and/or from differences between such interpolations under different 

3 0 physiological conditions or between normal and diseased tissue. 

The blood-nulling magnetic resonance sequence of FIGURE 2 is optionally 
repeated with a repetition time (TR) to acquire reconstructed images for multiple slices or for 
multiple three-dimensional (3D), two-dimensional (2D), or one-dimensional (ID) spatial 
encodings, for example to obtain a three-dimensional reconstructed image volume. In 
35 FIGURE 2, an exemplary succeeding spatially non-selective inversion recovery pulse 74 1 
corresponds to the beginning of the next repetition of the blood-nulling imaging sequence. 
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Optionally, because the blood nulling establishes a substantially negligible blood signal, it is 
possible to acquire multiple slices or multiple ID, 2D, or 3D spatial encodings per single 
nulling condition. This applies to both imaging and spectroscopy applications. 

Reconstructed images acquired using inversion-recovery blood-nulling 
5 magnetic resonance sequences such as the exemplary sequence of FIGURE 2 can be used in 
various ways. For example, while it is known in the art to process a reconstructed image to 
determine proton density (p), Tl, and T2 weightings, the determined weighting values are 
typically composite values having contributions from both blood and tissue. By applying 
substantially the same processing to the blood-nulled reconstructed images, more accurate p, 

10 Tl, and T2 weightings for tissue can be obtained without interference from the blood signal. 
Optionally, other tissues are nulled in addition to blood, to minimize contributions from those 
other tissues or components. For example, such combined blood and tissue nulling can be used 
to largely isolate a magnetic resonance signal from cerebral spinal fluid (CSF) in the brain. 

Reconstructed images acquired using inversion-recovery blood-nulling 

15 magnetic resonance sequences contain information pertaining to vascular space occupancy 
insofar as the images includes contributions from tissue but substantially exclude contributions 
from the blood volume. By applying a physiological perturbation, the blood volume can be 
changed. For example, cerebral blood volume undergoes vasodilation responsive to visual 
stimulation and breath-hold. Similarly, cerebral blood volume undergoes vasoconstriction 

2 0 responsive to hyperventilation. Blood volume changes can also be induced by administration of 
a selected drug, for example, but not limited to, for stress testing or the assessment of vascular 
compliance. Moreover, certain diseases such as cardiac ischemia, stroke, cancer, vascular 
deformations, and the like, represent a chronic or transient physiological perturbation that can 
cause a change in blood volume detectable with the methodology described herein. 

25 In most of the above cases, the change in blood volume occurs principally in 

the microvessels rather than in the large vessels and parenchyma. Advantageously, the blood 
volume changes measured using the blood-nulled reconstructed images reflect the parenchymal 
blood volume (denoted BV herein) which substantially corresponds to the volume of the 
microvessels without contributions from the larger vessels. In contrast, vascular volume effects 

30 measured by BOLD and many other existing techniques include the large vessels and 
parenchyma and other tissues close to these vessels (for example, cerebral spinal fluid). For 
parenchymal tissue (tissue with perfused blood), the parenchymal vascular space occupancy 
(VASO, also denoted £ herein) is given by: 
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BV + V. V 

13 v T v tissue v par 



where BV is the blood volume, V tissuc is the pure tissue volume (without blood), and V par is the 
volume of the parenchymal tissue, that is, Vpa^BV+V^sue. The measured parenchymal vascular 
5 space occupancy £ advantageously is more sensitive to physiological perturbation, including 
permanent disease-induced perturbation, than is the total vascular volume. Changes in 
large-vessel volume (that is, outside parenchymal regions) are also accessible by the blood 
nulling approach. At appropriate resolution, changes in large-vessel volume do not interfere 
with the indicated parenchymal blood volume changes due to the applied spatial encoding. 
1 0 For parenchymal tissue, the magnetic resonance signal (denoted by S herein) is 

proportional to a sum of the magnetization contributions of the microvessels and of the pure 
tissue. In general, the total signal S par from the parenchymal tissue (including both tissue and 
perfused blood) can be written as: 

1 5 S par = S bIood + S tissue (2), 



20 



30 



where S b i o0 d and Sti SSue are the blood and tissue signals, respectively. Inserting the vascular space 
occupancy £, and water density factors Cp ar and C b i 00 d for parenchymal tissue and blood, 
respectively, into Equation (2) yields: 

= C bIood •^M blood -e- TE ^- + (C par -^•C blood )-M tissue . e -™*- (3), 



where T2 b i 00 d and T2ti ssue are suitable time constants for the exemplary spin echo, readout 72 of 
FIGURE 2 due to the spin refocusing produced by the radio frequency pulse 84. For a gradient 
2 5 echo readout, the T2 values should be replaced by T2* values in Equation (3). The terms M b i 00 d 
and M tisS ue correspond to the initial transverse magnetization of the blood and the tissue, 
respectively, produced by the excitation radio frequency pulse 80. As a consequence of using 
the blood-nulling excitation sequence 70, however, the term M b i 00 d corresponding to the blood 
is substantially zero, and so Equation (3) reduces to: 



=(0^ -S-C^yM^ .e-™ 13 - (4). 
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A fractional parenchymal signal change AS/S occurring responsive to a physiological 
perturbation, defined as: 



5 



can be written by combining Equations (4) and (5) as: 

AS _ (Cp.-r ■C MBOd )-o- wna - -(C^ .^.e"™ 1 *- (6) 
S (C par -r st -C bIood )-e™^ 

10 where the superscript "act" references the perturbed state and the superscript "rest" references 
the unperturbed rest state or another suitable reference state. Dividing the numerator and 
denominator of Equation (6) by C par e*' re/T2(tisSttC ' rest) , the fractional parenchymal signal change 
AS/S can be written as: 



(1 _ £ act . ^blood ^ t e -TE/(l/T2fS uc -l/T2[S JC ) ^rest m C blood ^ 

15 = S= (7), 

S ^ j ^ rest ^ ^ blood ^ 

Cpar 

Making use of the series expansion approximation e*~l+x and defining an apparent change in 
transverse relaxation rate of tissue AR2 as: 



20 AR2 = ( 8 )> 

T9 act ^rest 

tissue tissue 



Equation (7) can be rewritten as: 



-AR2TE - ^ act > Cbl00d +^ act .£Mool. AR2 . T E + ^ rest .^blood 
AS C par C par C P ar ^ 

S (l-£ Test ^^l° od ) 



25 
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Defining a rate of change A£ of the vascular space occupancy as A£=4 act -4 rest and neglecting the 
term 2, act (C b iood/Cp a r) AR2 TE in Equation (9) and simplifying yields: 



-AR2-TE - A§ - — MqocL 



AS c par 

S (l-£ rest ) 
Cpar 



(10). 



FIGURE 3 shows a block diagram of an exemplary processor 100 that 
computes a rest blood volume BV"*' and a blood volume change ABV based on Equation (10). 
The processor 100 receives imaging data input from the image memory 46. The imaging data 
includes images taken at a plurality of echo times TE 102 with the imaging subject at rest, and 

1 0 corresponding images taken at the plurality of (that is, two or more) echo times TE with the 
imaging subject in a perturbed state. A processor 104 computes the fractional signal difference 
AS/S between the image of the subject in the rest state and the image of the subject in the 
perturbed state for each echo time TE to form a table 106 of fractional signal difference AS/S 
versus echo time TE. Recognizing that Equation (10) is a linear equation in echo time TE, a 

1 5 linear regression processor 110 computes a slope 112 and an ordinate-intercept 114 of the AS/S 
v. TE relationship. 

A rest blood volume BV** 1 calculator 116 computes the rest blood volume 
B yrest llg from the slope 112 according to: 



on 1 -AR2 -AR2 nn 

20 slope = - = — (11), 

i^rest .Cblood BV rest C blood 



C 



P ar v par c par 



where Equation (11) is the slope component of the linear relationship between AS/S and TE of 
Equation (10). The rightmost side of Equation (11) is written with the rest parenchymal 
vascular space occupancy £ rest replaced by blood volume BV rest divided by the volume V par of 
25 the parenchymal tissue under consideration in accordance with Equation (1). The BV rest 
calculator 116 uses literature values for the water density factors C par and C b i 00 d> the known 
volume V par of the parenchymal tissue under consideration, and a value for AR2 calculated 
from T2 dssue>BCt and T2 tissue , re st values obtained by fitting the signal value S at several echo times 
TE for each of the perturbed and rest states, respectively. If V par is unknown, the method can 

30 be used to determine the absolute blood volume fraction or vascular space occupancy (g rest ). 
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A blood volume change ABV calculator 122 computes the blood volume rate 
of change ABV 124 from the ordinate-intercept 114 according to: 



25 



Intercept = 



( ABV 


c ^ 

^ blood 


lv par 


C par J 


gyrest 


c ^ 

^blood 


^par 


C par j 



(12), 



where Equation (12) is the ordinate-intercept component of the linear relationship between 
AS/S and TE of Equation (10) written with the rest parenchymal vascular space occupancy 
replaced by BV rest /V par and the change in vascular space occupancy A£ is replaced by 
ABV/Vpar. The ABV calculator 122 uses literature values for the water density factors Cp ar and 

1 0 Cbiood, the known volume V par of the parenchymal tissue under consideration, and the rest blood 
volume BV^ 118 computed by the BV rest calculator 116. Alternatively, the absolute blood 
volume fraction or vascular space occupancy (<t rest ) can be used if V par is unknown. 

FIGURE 4 shows a block diagram of another processor 130 that computes an 
absolute blood volume BV. Computation of BV by the processor 130 is not based on a 

1 5 transient response to a physiological perturbation. Rather, data from the image memory 46, 
which may be a single image, is processed using Equation (4). The image is preferably 
acquired using a short echo time TE, so that the term exp(-TE/T2ti SSue ) in Equation (4) is 
suitably approximated as unity so that the signal S is given by: 

BV 

2 0 S par s (C par - — . C bl00d ) . M tissue (13), 

where the parenchymal vascular space occupancy % is replaced by the expression in Equation 
(1). The value of Mtissue at the inversion time TI 60 is related to longitudinal tissue relaxation 
time Tl tissue according to: 



Mtissue (TI) = l-2e- TI/Tltissue +e TR/T1 tissue (14)j 



30 



where TR is the repetition time of the pulse sequence (indicated, for example, in FIGURE 2). 
A suitable value for Tl tissue can be measured using a known magnetic resonance technique such 
as saturation recovery pulse sequence and inversion recovery pulse sequence or, preferably, by 
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a series of inversion pulses so that blood remains substantially nulled, while tissue signal 
decays with the representative Tissue- In the case of brain imaging, it is known that the Tl tissU e 
values for white and gray brain matter differ. Hence, as shown in FIGURE 4 a Tl W hi te value 132 
and a Tlgray value 134 are preferably determined. Notice that S par also contains instrumental 

i factors, for instance related to coil sensitivity and excitation homogeneity, the combined 
contribution of which we will call IF. This constant can also be determined during the Tl 
measurement from the voxels in ventricles, which only contains CSF and has a proton density 
very close to pure water. An Mfesue processor 136 computes corresponding Mdssue values at the 
inversion time TI 60 for white and gray brain matter, and a normalization processor 140 

) normalizes the image stored in the image memory 46 to produce a normalized image 142 
according to: 



S m (tone)- (15), 
IF *M 



5 where in Equation (15) the tissue magnetization is selected as the magnetization term 

corresponding to the Tl^y value 134 for gray matter, and is selected as the magnetization term 
corresponding to the Tl white value 136 for white matter. A blood volume processor 144 
estimates the absolute blood volume BV 146 or the absolute blood volume fraction or vascular 

space occupancy (% rest ) by applying Equation (13) to the normalized image 142. For pure 
0 parenchyma, an additional approach to determine blood volume is to use the tissue Tl value to 
null the tissue and only look only at remaining blood. Using the same IF determined in the 
blood nulling experiment, the signal in such a tissue-nulled voxel at short TE is given by: 



■••"•-iFfer 0 "--* (16) 



5 



from which the BV (Equation (1)) and/or vascular space occupancy can be determined for 
known Cbi 0 od< 

FIGURE 5 shows a processor 150 that computes a clinical image from 
blood-nulled magnetic resonance images. The processor 150 is suitable in diagnostic 
! 0 applications for detecting an abnormality such as vasodilation or vasoconstriction caused by 
certain diseases such as cardiac ischemia, stroke, cancer, and the like. The image memory 46 
stores an image of the imaging subject to be diagnosed. The image is compared with a 
reference image 152 by a difference processor 154 to generate a clinical image 156. 
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In one suitable approach, the difference processor 154 computes an absolute 
difference between the images, so that normal areas appear dark in the clinical image 156 
whereas abnormal areas appear brighter in the clinical image 156 (or vice versa) due to 
differences between the subject image and the reference image in the vicinity of the 
5 abnormality. Optionally, the difference processor 154 computes a signed difference with a 
constant intensity level offset. In this approach, regions of vasodilation and regions of 
vasocontraction have opposite intensity polarities respective to the constant intensity level. 

The reference image 152 can be obtained from various sources. For example, 
an image of a normal subject can be employed for the comparison. Alternatively, in the case of 

10 a bilaterally symmetric anatomical structure such as the brain, a contralateral image can be 
used. For example, a suspect right-side of the brain can be compared with a presumed normal 
left-side of the brain, preferably after a suitable left-right transposition of the contralateral 
comparison image. Still further, in cases of anatomical structures that are substantially 
homogeneous in the absence of an abnormality, a suspect portion of the organ can be compared 

1 5 with presumed normal portion of the same organ. 

In the case of brain diagnoses, the normalized image 142 of FIGURE 4 is 
optionally employed in the processor 150. The normalized image 142 advantageously has 
suppressed contrast between white matter regions and gray matter regions of the brain image, 
which makes the vascular space occupancy contrast more visible. 

20 In some clinical applications, large vessels changes may occur in addition to 

microvascular changes. These can also be detected, and, at sufficient spatial resolution, can be 
separated from the microvascular changes. 

The invention has been described with reference to the preferred embodiments. 
Obviously, modifications and alterations will occur to others upon reading and understanding 

2 5 the preceding detailed description. It is intended that the invention be construed as including 
all such modifications and alterations insofar as they come within the scope of the appended 
claims or the equivalents thereof. 
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Claims 

Having described the preferred embodiments, the invention is now claimed to be: 

1. A magnetic resonance method including: 

performing a blood-nulling magnetic resonance excitation sequence (70) that 
substantially nulls a magnetic resonance signal from blood; and 

subsequent to the performing of the blood-nulling magnetic resonance excitation 
sequence (70), performing a readout magnetic resonance sequence (72) to acquire a magnetic 
resonance signal from tissue other than the nulled blood. 

2. The magnetic resonance method as set forth in claim 1, wherein the performing of a 
blood-nulling magnetic resonance sequence (70) includes: 

performing an inversion recovery magnetic resonance excitation sequence (70) having 
an inversion time (60) to substantially null the magnetic resonance signal from blood. 

3. The magnetic resonance method as set forth in claim 2, wherein the performing of 
an inversion recovery magnetic resonance sequence (70) includes: 

applying an inversion radio frequency pulse (74); 
delaying for the inversion time (60); and 
applying an excitation radio frequency pulse (80). 

4. The magnetic resonance method as set forth in claim 3, wherein: 

the applying of the inversion radio frequency pulse (74) is performed without an 
accompanying spatially selective magnetic gradient pulse; and 

the applying of the excitation radio frequency pulse (80) is performed with an 
accompanying spatially selective magnetic field gradient pulse (82). 

5. The magnetic resonance method as set forth in claim 3, wherein the performing of 
an inversion recovery magnetic resonance sequence (70) further includes: 

applying additional inversion radio frequency pulses to maintain blood in a 
substantially nulled condition. 

6. The magnetic resonance method as set forth in claim 3, wherein the inversion radio 
frequency pulse (74) is a 180° pulse and the excitation radio frequency pulse (80) is a 90° pulse. 
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7. The magnetic resonance method as set forth in claim 2, further including 
determining the inversion time (60) for nulling blood by: 
determining a Tl value of blood; and 

selecting the inversion time (60) for substantial blood nulling based on the determined 
Tl value of blood. 



8. The magnetic resonance method as set forth in claim 7, wherein the determining of 
a Tl value of blood includes: 

acquiring a representative blood sample; and 

measuring the Tl value of the representative blood sample. 

9. The magnetic resonance method as set forth in claim 7, wherein the determining of 
a Tl value of blood further includes: 

determining a hematocrit of the blood; and 

correcting the Tl value of blood for the determined hematocrit 

10. The magnetic resonance method as set forth in claim 2, further including: 
identifying at least a main magnetic field strength and a repeat time; and 
determining the inversion time (60) for nulling blood based on a predetermined 

relationship (64) between the inversion time for nulling the blood and the identified main 
magnetic field strength and repeat time. 

11. The magnetic resonance method as set forth in claim 2, further including: 
optimizing an inversion time of a calibration inversion recovery magnetic resonance 

excitation sequence to minimize a magnetic resonance signal of a large blood vessel; and 

selecting (66) the inversion time (60) for nulling blood as the optimized inversion time 
of the calibration inversion recovery magnetic resonance excitation sequence. 

12. The magnetic resonance method as set forth in claim 2, further including: 
generating a reconstructed image from the acquired magnetic resonance signal. 

13. The magnetic resonance method as set forth in claim 12, further including: 
subsequent to performing the readout magnetic resonance sequence, inducing a 

physiological perturbation; 

subsequent to inducing the physiological perturbation, repeating performing the 
inversion recovery magnetic resonance excitation sequence (70) and performing the readout 
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magnetic resonance sequence (72) to acquire a second magnetic resonance signal from tissue 
other than the nulled blood; and 

generating a perturbation reconstructed image from the acquired second magnetic 
resonance signal. 

14. The magnetic resonance method as set forth in claim 12, further including: 
subsequent to performing the readout magnetic resonance sequence (72), inducing a 

physiological perturbation; 

subsequent to inducing the physiological perturbation, repeating performing the 
inversion recovery magnetic resonance excitation sequence (70); 

subsequent to repeating the performing of the inversion recovery magnetic resonance 
excitation sequence (70), performing a plurality of readout magnetic resonance sequences (72) 
each having a different echo time to acquire a plurality of magnetic resonance signals 
corresponding to the plurality of echo times; 

generating a plurality of perturbation reconstructed images from the acquired plurality 
of magnetic resonance signals corresponding to the plurality of echo times; and 

determining a temporal evolution of a physiological response to the physiological 
perturbation based on the plurality of perturbation reconstructed images. 

15. The magnetic resonance method as set forth in claim 14, wherein the determining 
of a temporal evolution includes: 

computing a change in vascular space occupancy signal between perturbation 
reconstructed images and corresponding unperturbed reconstructed images for each echo time 
to produce change in vascular space occupancy signal versus echo time data (106); and 

fitting the change in vascular space occupancy signal versus echo time data (106) to a 
mathematical model to obtain a blood volume parameter value (118). 

16. The magnetic resonance method as set forth in claim 12, further including: 
determining a Tl value of tissue (132, 134); 

computing a tissue magnetization (M) based on the Tl value of the tissue (132, 134); 

generating a normalized reconstructed image (142) by dividing the reconstructed image 
by the tissue magnetization; and 

estimating a blood volume parameter value (146) based on the normalized 
reconstructed image (142). 
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17. The magnetic resonance method as set forth in claim 16, wherein the reconstructed 
image corresponds to a subject brain region, and determining the Tl value of the tissue (132, 
134) includes: 

determining a first Tl value (132) corresponding to white brain matter, and 
determining a second Tl value (134) corresponding to gray brain matter; 
wherein the dividing of the reconstructed image by the tissue magnetization uses for 
each image element a selected one of the first and second Tl values (132, 134) based on a 
classification of local tissue type as one of gray brain matter and white brain matter. 

18. The magnetic resonance method as set forth in claim 12, wherein performing the 
readout magnetic resonance sequence (72) effects imaging of a subject brain region of a subject 
brain, and the method further includes: 

providing a reference image (152) of a reference brain region; and 
comparing the reconstructed image with the reference image (152) to detect an 
abnormality of the subject brain region. 

19. The magnetic resonance method as set forth in claim 18, wherein providing the 
reference image (152) of a reference brain region includes repeating: 

performing the inversion recovery magnetic resonance excitation sequence (70), 
performing the readout magnetic resonance sequence (72), and 

generating the reconstructed image on the reference brain region to generate the 
reference image (152). 

20. The magnetic resonance method as set forth in claim 19, wherein the reference 
brain region is selected from a group consisting of: 

a brain region of a contralateral side of the subject brain corresponding 
to the subject brain region, 

a brain region of a brain other than the subject brain which 
corresponds to the subject brain region, and 

a brain region of the subject brain other than the subject brain region. 

21. The magnetic resonance method as set forth in claim 12, further including: 
computing at least one of a proton density weighting of tissue, a Tl weighting of tissue, 

and a T2 weighting of tissue without interference from a magnetic resonance signal of blood 
based on the reconstructed image. 
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22. The magnetic resonance method as set forth in claim 2, wherein performing the 
readout magnetic resonance sequence (72) includes performing one or more of: 

a single-shot imaging sequence, 
a single-shot echo planar sequence, 
a multi-shot imaging sequence, 
a spectroscopy sequence, 
a multiple slice image, 

a one-dimensional, two-dimensional, or three dimensional spatial encoding sequence, 
a fractional k-space acquisition sequence, 
a spin echo readout sequence, and 
a gradient echo readout sequence. 

23. A magnetic resonance system including: 

a blood nulling means (10, 12, 16, 18, 30, 32, 34) for performing a blood nulling 
magnetic resonance excitation sequence (70) that substantially nulls a magnetic resonance 
signal from blood; and 

a readout means (10, 12, 16, 18, 30, 32, 36) for performing a readout magnetic 
resonance sequence (72) to acquire a magnetic resonance signal from tissue other than the 
nulled blood, the readout means (10, 12, 16, 18, 30, 32, 36) operating subsequent to operation 
of the blood nulling means (10, 12, 16, 18, 30, 32, 34). 

24. The magnetic resonance system as set forth in claim 23, wherein the blood nulling 
means (10, 12, 16, 18, 30, 32, 34) includes: 

an inversion recovery means (10, 12, 16, 18, 30, 32, 34) for performing an inversion 
recovery magnetic resonance excitation sequence (70) having an inversion time (60) in which 
magnetic resonance of blood is substantially nulled. 

25. The magnetic resonance system as set forth in claim 24, further including: 

a means (64) for determining the inversion time (60) based on a set of values including 
at least a magnetic field strength value and a repeat time value. 

26. The magnetic resonance system as set forth in claim 24, further including: 

a means (62) for measuring a Tl value of blood, the inversion time (60) being obtained 
from the measured Tl value of blood. 



27. The magnetic resonance system as set forth in claim 24, further including: 
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a reconstruction means (44) for generating a reconstructed image from the acquired 
magnetic resonance signal. 

28. The magnetic resonance system as set forth in claim 27, further including: 

a means (100, 130) for computing a blood volume parameter value (118, 146) from the 
reconstructed image. 

29. The magnetic resonance system as set forth in claim 28, wherein the means (130) 
for computing a blood volume parameter value (118, 146) includes: 

a means (136, 140) for normalizing the reconstructed image based on a Tl value of 
tissue (132, 134) to generate a tissue-normalized reconstructed image; and 

a means (144) for computing the blood volume from the tissue-normalized 
reconstructed image. 

30. The magnetic resonance system as set forth in claim 28, wherein the means (100) 
for computing a blood volume parameter value (118, 146) includes: 

a means (104) for computing an intermediate parameter (106) functionally related to 
blood volume for a plurality of reconstructed images produced by repetitively invoking the 
readout means and the reconstruction means with a corresponding plurality of echo times 
(102); and 

a means (110) for fitting a parameterized model to the intermediate parameters (106) 
and the corresponding echo times (102), the parameterized model having parameters including 
a rest blood volume (118) and a blood volume change (124). 



31. The magnetic resonance system as set forth in claim 27, further including: 
a means (154) for combining the reconstructed image with a reference image (152) to 
identify an abnormality in the reconstructed image. 
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